Interacting Geodesies: Binary Systems around a Black Hole. 
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We present a novel method to study interacting orbits in a fixed mean gravitational field associated 
with a solution of the Einstein field equations. The idea is to consider the Newton gravity among 
the orbiting particles in a geometry given by the main source. We apply the technique in the of 
study two and three self-gravitating particles moving around a black hole, i.e., in a Schwarzschild 
geometry. We also compare with the equivalent Newtonian problem and noted differences in the 
structural stability, e.g., binary systems were found only in the general relativistic approach. 
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General relativistic effects are usually neglected in astrophysical systems in which many gravitational sources are 
,-0 ■ considered. In most of the known examples these effects are insignificant. Also the n-body problem is intractable in 
Einstein theory without large approximations. Nevertheless the difference between the Newtonian and the General 
Relativistic approaches becomes important close to compact objects and around very massive structures. When non 
interacting particles are been considered, the problem can be solved computing time-like geodesies. When more than 
one source is considered (in the non-stationary case) numerical techniques might be used to solve the Einstein field 
I \ equations. However finding numerical solution to the simple two body problem is already a non trivial problem, see 
for instance [lj. In some cases, analytical and numerical solutions to the n-body problem in GR are found as sequences 
of post-Newtonian approximations ■ Following the post-Newtonian idea, there are some predictions of relativistic 
effects using a relativistic celestial mechanic Other approximation techniques are also available, the simplest 
consists in solving a typical Newtonian problem using pseudo-potentials that simulates the GR aspects relevant in 

■ the studied situation UQ- 

^C) \ We suggest in this article a novel technique that makes possible the study of small massive objects orbiting around 

■ a fixed structure with a large mass. The idea is to consider the Newtonian gravity among particles evolving in a 
**^h geometry associated with a solution to the Einstein field equations. 

We construct within the Newtonian gravity a quadri-force. Then this force is used to modify the geodesic equation 
to allow a "Newtonian" interaction among them. In this scenario the damping due to the emission of gravitational 
radiation is not considered. 

In the context of Schwarzschild geometry, we study planar interacting orbits and compare the results with the 
equivalent situation in the usual Newtonian theory. We see a great difference in the structure stability of the studied 
systems. The most relevant result is that in our approach the generation of binary systems is favoured if compared 
with the equivalent pure Newtonian situation. In n-body simulations the production of binary systems is a key issue 
in the formation of astrophysical structures as galaxies and star clusters |fjj . 

In our proposition the interacting particles do not change the mean gravitational field. In the general relativistic 
context this field is represented by a solution to the Einstein equations. Non- interacting particles move following 
geodesies, i.e., their world lines obey the geodesic equations, 

i A V A i^ = + T" ap x a xP = 0, (1) 

where = g^ 1 (g a -yj3 + gpj.a — 9a/3,j), and g a p is the metric tensor for the known solution. The dots represent 
ordinary derivation with respect to the proper time r and commas denote the usual partial derivation with respect 
to the coordinate. 

The next step is to use Newton gravity to compute the interaction among orbiting particles. The Newtonian 
gravitational potential generated by a set of N massive particles is 

N n 

^|r-r a | 

In the rest of the article we shall use geometrical units: the Newton constant G = 1, the light speed c = 1, we also 
assume that the black hole mass is unit (M = 1). 

Therefore the force f that acts on the particle b due to the other orbiting particles is 

f (6) = -m 6 V 6 $(r a f r b ). (3) 



I 

o 

H 



2 



Now we look for a quadri-force that represents in some way the Newtonian interaction among the particles in a fixed 
pseudo-Riemannian geometry. The spatial components of the quadri-force is obtained from f given in J^J. From 
f = dp/dt and F l = dp l /dr we get F l = uq/ 1 , where u 11 = dx^/dr. We want a quadri-force F^ orthogonal to the 
quadri- velocity u M . Hence, we use u^F^ = uqF° + uoUif 1 = to obtain F° — —Uif 1 . We can write F M in a more 
adequate form as 



The signature of guv is taken as —2. 

Finally, we introduce F^ in the right side of the geodesic equations representing an external force to define the 
motion equations for a particle of mass m a , 



From now on we consider that all the particles moving around the center of attraction have the same mass m. 

First we consider the limits of (jHJ) for large and small masses. When the mass of the interacting particles m is much 
smaller than the fixed source mass M, we have nothing but the geodesic equations for a set of test-particles. For 
ra » M we get the special relativistic analogous to the Newton second law. And for small speeds, i.e., uq — > 1 and 
Ui — » 0, the Newtonian motion equations for a system of iV gravitating particles are recovered. 

The method is "self-inconsistent" since some of the mathematical principles of a general relativistic framework are 
not satisfied. For instance, we are using an instantaneous interaction among the particles. The distance between 
two particles calculated in the Black Hole framework was used to compute the force of attraction without correction 
related to the different speed of the particles. The force is not covariant. 

We must have these imperfections in mind to determine the limitation of the method and to check numerical results. 
Since the geometry will be kept fixed the mass of the orbiting particles must be much smaller than the main source. 
We have to systematically control variations of the total angular momentum and energy of the orbiting particles, they 
should be constant. 

We use the above described approach to study particles moving around a spherical symmetric source, i.e, we 
numerically solved © when the field g^ u is the Schwarzschild metric. The initial conditions are chosen to have for 
each particle a timelike worldline and that the spatial projection of these worldlines represent a bounded system of 
particles moving on the plane $ = 7r/2. 

In the first examples, three particles moving close to the black hole are studied. For this purpose we take : L\ — 
6.752, E 1 = 0.988; L 2 = 6.937, E 2 = 0.989 and L 3 = 5.895, E 3 = 0.984. [In the Schwarzchild metric: L = r 2 0sin 2 tf 
and E = t/(l — 2M/r)]. Their initial distances from the center of attraction are, respectivily, r% — 42.5, r 2 = 45.0 and 
r 3 = 31.5 in geometric units. 

When we compute the orbits for three particles interacting only with a center of attraction {Fm = F^) = ^(3) = ^) 
this is equivalent to chose m = 0, as above mentioned we take the black hole mass as one (M = 1). For the already 
presented values of L and E we have bounded geodesic motion around the black hole FigQ] In FigElwe present three 
particles with equal masses m = 5 X 10~ 4 . The trajectories are quite different from the geodesies but the particles still 
have bounded motion. In Fig|3]we present the motion for a larger value of the masses, m = 3 x 10 -3 . In this case, 
one particle escapes (the dashed orbit) while the other couple form a binary system that falls into the black hole. 

Now we consider only two particles farther from the black hole than the preceding ones. For this purpose, we 
consider L% — 28.10, E\ = 0.99937 and L 2 = 28.16, E 2 = 0.99937. As in the previous examples, the non interacting 
orbiting particles follow bounded geodesic motions, Fig0] In FigEJwe present the trajectories for the pair of orbiting 
particles with m = 10~ 5 . In this case, we can clearly dintinguish the motion of each particles. In the second case 
is considered m = 10~ 4 , now the particles move in the same region of the space in such way that their orbits fill a 
thick ring, FigEl Finally, in Fig[7| the particles mass are m = 5 x 10" 4 and we have the formation of a stable binary 
system that remains in a confined motion around the black hole. 

It is instructive to compare our results with the equivalent pure Newtonian three-body problem. The gravitational 
potential field of the system in this case is 



mu g 




(4) 
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With mi = m 2 = m the motion equations of the particles are, 




(7) 
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These equations are numerically solved by using parameters that confine the particles in a motion around the central 
mass. The Newtonian energy per unit of mass, EN ewt , is related to the relativistic specific energy by E = sfl — 2E^ ewt . 
The specific angular momentum L and energy E used in the simulations are L\ — 6.55, E\ = 0.988; Li = 6.75, E2 = 
0.989, and L3 = 5.68, E3 = 0.984. These values are equivalent to the ones used in the first general relativistic 
simulations (Figs^ |21and[3J). 

We start with the simulation of three test-particles (no interaction among them) moving around an attraction center 
in elliptical orbits, FiglHl We present in FigOU interacting particles with m = 5 x 10~ 4 . We see that one of them 
quickly unbound and the remaining couple stay in a bounded motion for some time. With a stronger interaction, 
to = 3 x 10~ 3 presented in Fig llOl two particles quickly move far from the center of attaction while the remaining one 
stabilizes in a bounded elliptical motion closer to the black hole. 

Now we present results of two particles moving farther from the center of attraction in a situation similar to 
the analogous relativistic presented in FigQ] [5] and [7] For this aim, we use L\ = 28.10, E\ = 0.99937 and 
L-2 = 28.12, E2 — 0.99937. The motion of test particles is presented in Fig. El In FigElwe have m = 10 -4 the 
particles stabilize in quasi-elliptical orbits. One closer and other more distant to the center of mass if we compare 
with the non-interacting situation. For a stronger interaction, m = 3x 10~ 3 , we have a interesting situation as shown 
in Fig ll3l The particles alternate quasi-elliptical motions closer an farther to the center of attraction, i.e., while one 
of them are orbtitating in a large ellipse the other is in a small one. We make some simulations considering that the 
central mass is not fixed and has the unit mass (M = 1). No significative difference was noticed. 

Before analysing the results, it is important to explain some numerical cautions necessary in these general relativistic 
simulations. In the first examples (Fig|2and FigEJ the particles were very close to the black hole (remember that if 
the BH have 10M© the distance of 50 in geometric units will be about 50km) therefore the system is more unstable, 
see for instance 0, , and the error due to the inexact approach is larger since the general relativistic effects are very 
important. We use, to control such problems, the numerical values of the total angular momentum and energy of 
the orbiting particles (that should be constant) and their physical time - if they were very different we would have 
problems in the syncronization of the Newtonian interaction. The numerical variation of the total energy was always 
smaller than 1%, the total angular momentum varied at most 10% and the time less than 5%. We consider that it is 
the worst situation that we still can apply this technique, in the other examples with the general relativistic approach, 
FigEJand Fig. [HI trie deviations of the total angular momentum, total energy and time were very small - always less 
than 1%. 

By comparing the orbits obtained in our GR framework with the equivalent situation in Newtonian gravity we 
conclude that the strutuctual stability of the same system changes substantially with different approaches. The most 
interesting difference is the formation of binary systems noticed in the Fig|3| (a binary system that is captured by the 
black hole) and in the Fig0 (the binary bounded motion stabilizes around the black hole) . An analogous situation is 
very hard to be obtained in the Newtonian approach so close to the center of attraction. When the orbiting particles 
are very far from the main gravitational source, the Newtonian and the relativistic approach lead to similar results. In 
these cases, the gravitational force exerced by the central mass varies much more slowly than the interaction between 
the two orbiting particles. 
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FIG. 1: Geodesic motion of three particles around the black 
hole. The initial angular momentum and energy of each 
particle are, respectively, L\ — 6.752, E\ = 0.988; L 2 = 
6.937, E 2 = 0.989 and L s = 5.895, E s = 0.984. 
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FIG. 2: Now we use the same values of angular momentum 
and energy of the first figure but now with interacting parti- 
cles. They have the same mass, m = 5 x 10~ 4 . We obtain 
three particles in a bounded motion around the black hole. 
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FIG. 3: Using the same values of angular momentum and 
energy of the preceding figure and rn — 3x 10 -3 we have now 
that two of the particles form a binary system that falls into 
the black hole while the other stabilyzes in a bounded motion 
around the center. 



-500 500 

X 

FIG. 4: Geodesic motion around a black hole of two particles 
with Li = 28.10, E 1 = 0.99937 and L 2 = 28.16, E 2 = 0.99937. 
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FIG. 5: We maintain the values of angular momentum and 
energy of the preceding figure with two interacting particles 
with m — 10 . They have a bounded motion close to the 
geodesies. 
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FIG. 6: Using the same L and E of the last figure but with 
m = 10 -4 , we have that the motion of the particles is still 
bounded but now they are more irregular and we cannot easily 
distinguish them. 



■5 




FIG. 7: Increasing the mass to m — 5 x 10 4 we obtain a 
binary system moving around the center of attraction. 



FIG. 10: Increasing the mass, m = 3 x 10~ 3 , we have that 
two of the particles quickly escapes from the system while the 
remaining one tends to stabilizes in a elliptical motion. 




FIG. 8: Newtonian analogous to FigQ Three test-particles 
with Li = 6.55, Ei = 0.988; L 2 = 6.75, E 2 = 0.989. and L 3 = 
5.68, E 3 = 0.984. 




FIG. 11: The Newtonian analogous to the system presented 
in Fig^] is reproduced here. Two test-particles with Li = 
28.10, Ei = 0.99937 and L 2 = 28.12, E 2 = 0.99937. The 
Newtonian analogous to the system presented in Figgis re- 
produced. 
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FIG. 9: With the same values of L and E of the preceding 
figure, and rn — 5 x 10 -4 we get the Newtonian analogous to 
Fig|5| We can see that one of the particles escapes. 



FIG. 12: We use the same values of angular momentum and 
energy of the preceding figure with interacting particles with 
m — 10 -4 . The Newtonian analogous to the system presented 
in Figgis reproduced. We obtain, as a final configuration, 
two quasi-elliptical motion. 
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FIG. 13: Using the same values of L and E as in the previous 
figure and increasing to m — 1CF 3 we obtain that the motion 
of the particles are linked in such way that they alternate 
between a small and a large precessing ellipse. 



